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an external magnetic field similar to the Little-Parks effect and formation of mul- 
tiquantum vortex structures. The structure of a superconducting state is studied 
both analytically and numerically. 
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Little-Parks effect Q], i.e. oscillations of the critical temperature T c of multiply- 
connected superconducting samples in an applied magnetic field H, is one of the 
striking phenomenon demonstrating coherent nature of superconducting state. Such 
oscillatory behaviour of T C (H) is not a distinctive feature of a superconducting thin- 
wall cylinder and can be observed also in superconductors with columnar defects 
and holes (see references 01 H]) and mesoscopic simply-connected samples of 
the size of several coherence lengths [HI EJ- Generally, the oscillations of T c with 
a change in the external magnetic flux are caused by the transitions between the 
states with different vorticities (winding numbers) characterizing the circulation of 
the phase of the order parameter. For a system with cylindrical symmetry the 
vorticity parameter just coincides with the angular momentum of the Cooper pair 
wave function. The states with a certain angular momentum m can be considered as 
m-quanta vortices. Experimental and theoretical investigations of these exotic vortex 
structures (multiquanta vortices and vortex molecules) in mesoscopic superconductors 
have attracted a great deal of attention. As we change an external homogeneous 
magnetic field multiquanta vortices and vortex molecules can transform one into 
another via first or second order phase transitions. 

In this paper we focus on another possibility to create multiquantum vortex states: 
nucleation of superconducting order parameter in a hybrid system consisting of a thin 
superconducting film and an array of magnetic nanoparticlcs. The interest to such 
structures is stimulated by their large potential for applications (e.g., as switches 
or systems with a controlled artificial pinning). The enhancement of the depinning 
critical current density j c has been observed experimentally for superconducting 
films with arrays of submicron magnetic dots El E]> antidots 10 , and for 
superconductor-ferromagnet (S/F) bilayers with domain structure in ferromagnetic 
films 11 . The matching effects observed for magnetic and transport characteristics 
were explained in terms of commensurability between the flux lattice and the lattice 
of magnetic particles. Vortex structures and pinning in the S/F systems at rather 
low magnetic fields (in the London approximation) have been analysed in papers 

H3 ca ii3i csi nsi nzi na nni Eoi m] 

Provided the thickness of a superconducting film is rather small as compared 
with the coherence length, the critical temperature of the superconducting transition 
as well as the structure of superconducting nuclei should be determined by a two- 
dimensional distribution of a magnetic field component B z (x,y) (perpendicular to 
the superconducting film plane) induced by the ferromagnetic particles. Obviously, 
the highest critical temperature corresponds to the nuclei which appear near the 
lines of zeros of B z due to the mechanism analogous to the one responsible for the 
surface superconductivity (see, e.g., reference [22)) and domain wall superconductivity 
[23 El ESI- Provided these lines of zeros have the shape of closed loops, the winding 
number of a superconducting nucleus will be determined by the magnetic flux through 
the loop. Thus, changing this flux (e.g., increasing an external H field applied along 
the z axis) we can control the winding number. The resulting phase transitions 
between the multiquantum states with different m can cause the oscillations of T c . 
Such oscillatory behaviour has been, in fact, observed in reference j2E] for a Nb film 
with an array of GdCo particles. Note that a change in a slope of the phase transition 
curve T C (H) (which is probably a signature of the transition discussed above) has been 
also found in reference [23 for a Pb film with CoPd particles. Provided the dimensions 
of the sample in the (xy) plane are compared with the coherence length, we can expect 
a rather complicated picture which is influenced both by the sample edges and by the 
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distribution of an inhomogeneous magnetic field. For a model step-like profile of the 
magnetic field the resulting phase transitions between different types of exotic vortex 
states in a mesoscopic disc have been studied numerically in reference [25]. The 
interplay between the boundary effects and magnetic field inhomogeneity influences 
also the formation of multiquantum vortex states around a finite size magnetic dot 
embedded in a large area superconducting film [23 I3()j . The transitions between 
different multiquantum vortex states with a change in magnetic field and magnetic dot 
parameters were studied in references [2B1GII1IIIII f° r certain temperature values. These 
effects are closely related to the ones observed in mesoscopic and multiply-connected 
samples, and consequently, we can expect that the oscillations of T C (H) (analysed 
below) should also be a common feature of multiply-connected superconductors and 
thin-film systems with magnetic dots. 

In this work we do not consider the magnetic phase transitions in the mixed 
state for T < T c and focus on the oscillatory behaviour of T C (H) in a large 
area superconducting film caused only by the quantization associated with the 
characteristics of the inhomogeneous magnetic field produced by ferromagnetic 
particles. We neglect the influence of the edge and proximity effects in the S/F system 
and consider a nanoparticle only as a source of a small-scale magnetic field. Our further 
consideration is based on the linearized Ginzburg-Landau model: 

Here ^(r) is the order parameter, A(r) is the vector potential, B(r) = V x A(r), <I>o 
is the magnetic flux quantum, £(T) = £o/\/l ~~ T/T c0 is the coherence length, and 
T c o is the critical temperature of the bulk superconductor at B = 0. For the sake 
of simplicity we neglect the interference effects between the superconducting nuclei 
appearing near different nanoparticles (i.e. assume the interparticle distance to be 
rather large as compared with the superconducting nucleus size) and consider a single 
magnetic particle with a fixed magnetic moment chosen perpendicular to the film 
plane xy. For a rather thin film (of the thickness less than the coherence length) we 
can neglect the influence of the field components -B x , B y in the film plane and consider 
an axially symmetrical two-dimensional problem Q in the field B z (r) = H + b(r), 
where b(r) is the z-component of the field induced by the ferromagnetic particle and 
(r, 9, z) is a cylindrical coordinate system. Choosing the gauge Ag(r) = Hr/2 + a(r) 
one can find the solution of the equation |QJ in the form \&(r) = g m (r) exp(im#) / \/r , 
where m is the vorticity, and g m (r) cx r \ m \+ 1 /' 2 f or r _> o. The function g m (r) should 
be determined from the equation: 

" ($(r)/$ + m) 2 _ J_~ 



d 9rn 

dr 2 



9m — £2(rp^9 m ' W 

Here $(r) = 2nrAg(r) is the total flux through the circle of radius r. The lowest 
eigenvalue I /£ 2 (T) of the Schrodinger-like equation (J2J defines the critical temperature 
T c of the phase transition into a superconducting state. 

Obviously, for rather small fields H the superconducting order parameter can 
nucleate either far from the magnetic particle (r — > oo) where the critical temperature 
T C H is defined by the homogeneous field B z — H or in the region close to the circle 
of the radius ro where B z (ro) = and T c is controlled by the slope of B z (r) at 
r = ro and by the flux through the area of the radius ro. In the first case we obtain 
1 — T C H /T c0 = 27r|iJ|^Q/$ . For the second case we can analyse the behaviour of T C (H) 
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assuming that the characteristic length scale £ of the order parameter nucleus is much 
less than the characteristic scale of the magnetic field distribution. Within such local 
approximation (similar to the one used in reference |25j for the description of domain 
wall superconductivity) we can expand the flux in powers of the distance from rg: 



$(r) 
5>n 



irr Q B' z (r ) 
$ 



(r-r ) 2 



This local approximation is valid under the following conditions: 



B' z (ro) 



< 1 



and — < 1. 

ro 



Introducing a new coordinate t = (r — ro)/£ we obtain the dimensionless equation 

-^ + (t 2 -Qfg = Eg, 
where the parameters E and Q are given by the expressions: 



(3) 



E = e 

so 



1 



T 

T c o 



$ 



*\B' B (r )\ 



Q 




(4) 



We obtain E{Q) ~ Q 2 + % /^2Q when Q < -1, and E{Q) ~ when Q > 1. 
The minimal value of i£(Q) is E 1 = _E m i n — 0.904 at Q ~ 0.437. The final expression 
for the critical temperature reads: 



1 - 



c() 



J? 2 



min _B — 




(5) 



The superconducting nuclei are localized near the ferromagnetic particle at a distance 
ro. The states with different energetically favorable winding numbers m correspond to 
the multiquantum vortex structures very similar to the ones observed in a mesoscopic 
disc. As we change an external field H, we change the flux <&(ro) and, thus, change 
the energetically favorable vorticity number and position of the nucleus. 

To investigate the details of the oscillatory behaviour discussed above we consider 
a particular case of a small ferromagnetic particle which can be described as a point 
magnetic dipole with a magnetic moment M = Mzq placed at a height h over the 
superconducting film. The corresponding expressions for the field and the vector 
potential are: 



b(r) 



M(2h 2 



(r 2 + /i 2 ) 5 / 2 



a(r) 



Mr 



( r 2 + ^2)3/2 



(6) 



Introducing f m (f) = 9m{f)l \fr and a dimensionless coordinate p = r/h we obtain 
the equation J5Jl in the form: 

2 



1 d_ 

pdp 



dfm 



3V3 



P 



H 



(l+p 2 )3/ 2 




,(7) 



where N{ = AttM/ (3\/3h&o) is the dimensionless flux through the area with the 
positive field b(r) and bo — 6(0) = 2M/h 3 . In the limit of small fields H — > the 
nucleation of superconductivity occurs at large distances p and critical temperatures 
for different winding numbers m are very close. Thus, in this limit the critical 
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temperature is equal to T C H and is not sensitive to the presence of the dipole. 
Below this temperature we obtain a lattice of singly quantized vortices (with the 
concentration determined by H) which is surely disturbed under the dipole. Note, that 
the behaviour of T c in this low field regime should modify provided we take account 
of a finite distance between the magnetic particles. For large absolute values \H\ 
(much larger than the maximum field induced by the dipole) we obtain the following 
asymptotical behaviour of T c : 1 — T c /T c q = 2tt^q(—H — bo)/$o for negative H and 
1 - T c /T c o = 2<o(-ff - &o/(25V5))/$ for positive H values (here -6 /(25V5) is the 
minimum of the dipole field). The superconductivity nucleates near the minima of 
the total field \B Z \ and, thus, is localized near the dipole. In the intermediate field 
region (— 1 < H/bo < l/(25 v / 5)) we should expect the oscillatory behaviour of T c 
discussed above. The number of oscillations is controlled by the parameter N[. We 
have carried out the numerical calculations of equation (JZJ) for the various JVf values. 
For the numerical analysis of the localized states of equation Q we approximated 
it on a equidistant grid and obtained the eigenfunctions f m (p) and eigenvalues by 
the diagonalization method of the tridiagonal difference scheme. The results of these 
calculations as well as the analytical dependence of T c given by the expression JSJ are 
shown in figure ^ 
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Figure 1. Critical temperature as a function of an external magnetic field for 
Ni = 4 (a) and iVf = 10 (b). Solid line is a result of direct numerical simulations 
of equation J7J. Dash line is obtained from the analytical formula JFJ. Certain 
winding numbers m for different parts of the phase transition line are shown. 

We observed a remarkable asymmetry of the phase transition curve (T C (H) ^ 
T c (—H)) which is caused by the difference in the distributions of positive and negative 
parts of the dipole field b(r): the maximum positive field (bo) is much larger than 
the absolute value of the minimum negative field (6o/(25v / 5))- As a result, the T c 
oscillations appear to be most pronounced for negative H which compensate the 
positive part of the dipole field. Taking M ~ 3 X 10~ n G ■ cm? (for a ferromagnetic 
particle with dimensions 300 nm x 300 nm x 300 nm and magnetization ~ 10 3 G), 
h ~ 300 nm we obtain Nf ~ 10, bo ~ 10 3 G and the characteristic scales of T c 
oscillations AH - 100 Oe, AT C - 10~ 2 T c0 ~ 0.1 K for a Nb film with £ ~ 40 nm 
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and T c0 ~ 8 K. 

We expect that the oscillatory behaviour of T c can be observable, e.g., 
in magnetoresistance measurements of thin superconducting films with arrays of 
ferromagnetic particles. The superconducting nuclei localized near the particles should 
result in the partial decrease in the resistance below the oscillating T C (H). As we 
decrease the temperature below T C (H) the supeconducting order parameter around a 
single particle becomes a mixture of angular harmonics with different m values and 
we can expect the appearance the phase transitions similar to the ones discussed in 
papers |281 15U| . With the further decrease in temperature the whole film becomes 
superconducting and resistivity turns to zero. 
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